Robust self-testing quantum states and binary nonlocal XOR games 
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Self-testing a quantum state means verifying the existence of a certain multipartite quantum state 
based only on knowledge of the probability distributions of certain measurements on the state. Self- 
testing is possible even when the measurement devices are untrusted. This paper is concerned with 
proving robust (i.e., error-tolerant) self-testing results. We devise a general scheme for proving such 
results for nonlocal binary XOR games. We apply this scheme to prove a result for the Greenberger- 
Horne-Zeilinger paradox. 



Suppose that two parties, Alice and Bob, share a bipar- 
tite quantum state system {Qa,Qb) and each possesses 
a prescribed set of measurements on his or her subsys- 
tem. Suppose that the state of the system is not knovirn, 
but suppose that we know the joint probability distri- 
butions of Alice's and Bob's measurement-outcomes for 
any pair of prescribed measurements. In some cases, this 
information alone is sufficient to deduce the state of the 
system {A,B). If a system of two qubits is such that 
joint measurements along the x-axis will always agree, 
and joint measurements along the z-axis always agree, 
then the state of the system must be the Bell state 

Taking this idea a step further, we may ask whether 
we can determine the state of a quantum system based 
only the outcome distributions of some unknown quan- 
tum measurements. This is the notion of self-testing of 
quantum states, and it is possible in some cases. For ex- 
ample, if a collection of binary measurements on a system 
is known to achieve the optimal quantum violation of the 
CHSH inequality, then the system must be equivalent to 
a collection of $+-pairs 

The notion of self-testing of quantum states is cru- 
cial for quantum protocols that involve untrusted devices. 
For example, in the task of randomness expansion from 
an untrusted device, a party wishes to use an untrusted 
quantum device to generate random numbers. The party 
has no knowledge of the inner workings of the device, and 
can base his trust in the randomness of the device only 
on its observed classical outputs. Self-testing results for 
quantum states give the basis for security proofs for such 
devices. The idea of self-testing in the context of secu- 
rity proofs goes back at least as far as Mayers and Yao 
in 1998 and it can be seen in more recent work on 
randomness expansion 0, 

In order for self-testing results to be practically use- 
ful in randomness expansion protocols, they must tol- 
erate error. It is not sufficient to know that particular 
measurement-outcome distributions guarantee the exis- 
tence of a particular state we must also know that 
when the distributions are close to some prescribed dis- 
tributions, then the underlying state must be close to the 



target state 

McKague et al. proved that the CHSH paradox is 
a robust self-test for a Bell state [S]. In independent 
work, the authors of the present paper proved that the 
Greenberger-Horne-Zeilinger (GHZ) paradox is a robust 
selttest for the GHZ state ^ (|000) + |111)) 3. The 
purpose of the present paper is to provide a general 
scheme for proving robust self-testing results which ap- 
plies to both the GHZ and CHSH games. The scheme 
builds on the method used in [6| . 

We work within the class of nonlocal binary XOR 
games (that is, binary games for which the scoring rule 
is based only on the XOR of the outputs). We show 
that if a binary nonlocal XOR game satisfies certain sim- 
ple criteria-based on the behavior of an associated mul- 
tivariable sinusoidal function-then it satisfies a robust 
self-testing result. The criteria that we prove is fairly 
easy to use provided that the number of parties n is not 
too large. As a demonstration of the theory, we provide 
a short proof of a robust self-testing result for the GHZ 
paradox. 

We note that a robust self-testing result for the GHZ 
gradox was obtained independently in the recent paper 

Definitions. For our purposes, a binary nonlocal XOR 
game is simply a function /: {0, 1}" — )■ R. The function 
/ describes a scoring rule for the game: if the input se- 
quence is («i, 12, • ■ • ,in)i and the output sequence satis- 
fies ®feOfe = 0, then the score is /(ii, 12, . . . , i„); if the 
input sequence is (zi, 12, . . . , in) and the output sequence 
satisfies ®kOk — 1, then the score is —f {11,12, ■ ■ ■ ,in)- 

To any nonlocal game /, let us associate a polynomial 
Pf : C" — > C hke so: for any n-tuple (Ai, . . . , A„) of com- 
plex numbers, let P/ (Ai, . . . , A„) be equal to 

f{lu...,Zn)XY\^---K-. (1) 

(ji,...,i„)e{o,i}" 

For example, if g is the Clauser-Horne-Shimony-Holt 
(CHSH) game (5(1, 1) = -1, .g(0, 0) = g{0, 1) = 5(1, 0) - 
1) then 



P„ = 1 + Al + A2 - A1A2 



(2) 
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Additionally, for any binary nonlocal XOR game 
/ : {0, 1}" — )■ M, and any real numbers 9o,9i, . . . ,9n, let 
Zf{Oo, . . . , 9n) denote the quantity 



^ /(ii,...,«„)cos U'o + ^ife^fcj 

(ifc)e{o,i}" \ k J 

Thus, for example, if g is the CHSH game, then 
^3(^0,^1,^2) =cos(0o)+ cos(0o + 01 ) 



+ cos(0o 



cos^Wo 



2). 



(3) 



(4) 
(5) 



Note that the function Zf is 27r-periodic in every variable. 

The two quantities P/ and Zj are related by the fol- 
lowing identity. 

Z/(0o, . . ■,9„) = Re[e''oPf{e''\. . . , e'^")]. (6) 

Note also that 

\Pf{e''\...,e''-)\^ max 0i, . . . , (7) 

itz [— tt.tt] 

Quantum strategies. For our purposes, a quantum 
strategy for a binary n-player nonlocal game is a pure 
state 



e Qi ® Q2 ' 



(8) 



where each Qi is a Hilbert space, together with two pro- 
jective measurements 



(9) 



on the space Qj . These measurements can be more com- 
pactly expressed as Hermitian operators: 



M, 



(1) 



p(i-+) _ p(i -) 



(10) 
(11) 



J ■ 3 

The score for such a strategy is 

(V^l /(zi, . . . ,z„)Mf^^ M^) IV) . (12) 

Let us use the term nondegenerate n-qubit strategy to 
refer to a strategy whose Hilbert spaces Qj are all copies 
of and whose projection operators pj*'*-* are all one- 
dimensional projectors. 

Proposition 1. Let f : {0,1}" — > K 6e a nonlocal bi- 
nary XOR game. Let Q f denote the optimal score for f 
that can be achieved by a nondegenerate n-qubit .strategy. 
Then, 

Qf = ^^^ max |Py(Ai,...,A„)| (13) 

|Ai| = ... = |A„| = l 



Proof. Let ^^/i, llMj^-*, Afj^''}| J be a nondegenerate n- 

qubit strategy for /. Each of the operators Mj*-* is a 
Hermitian operator on a 2-dimensional space that has 
eigenvalues in the set {—1,-1-1}. After an appropriate 
change of basis, we may make the assumption that 



1 

1 



e'"^ 



(15) 



for some 0o, • • ■ , G [— tt, tt]. 

The score for this quantum strategy is clearly bounded 
by the operator norm of the operator 

M := ^ /(zi, . . . , tn)Mi''^ <E>...<E> M^^") (16) 

The operator M is on a Hilbert space which has basis 
{|aia2 . . . a„) | G {0, 1}}. If we take the elements of 
this basis in lexicographical order, the resulting matrix 
expression is a reverse-diagonal matrix: 







* 

* 



* 

* 







(17) 



The entries along the reverse diagonal are given by the 
expressions 



P 



(18) 



for (flfc) e {0, 1}". 

Using the simple observation that the eigenvalues of 
any matrix of the form 



zi 



Z2 



Z2 



Z\ 



(19) 



are ± j^ij ,± \zi\ , . . . ,± |z„|, we find that the operator 
norm of M is 



max 

(ai)e{0,l}'' 



P/ 



i(-i)""e. 



(20) 



Qi 



max Zf{UQ, 
9o,...,e„e[—n,7T] 



(14) 



Formula ([13)1 follows. Formula ([T4| follows also via equal- 
ity dS]). □ 
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Approximation. Let us say that two quantum strate- 
gies (V',{{Mf \Mj^)}},) and (0, {{ivf \ ivj^^}},) which 
operate on the same Hilbert spaces Qi,...,Qri are S- 
close if 



N. 



(0 



< S 



Vie{l,2,...,n},ie{0,l}, 



and 



Let us say that a nonlocal game / satisfies a second- 
degree approximation property if the following holds. 

• There exist constants C, eg such that any nondegen- 
erate n-qubit strategy whose score is within e < eg 
of the optimal score is (C-y/e)-close to an optimal 
nondegenerate n-qubit strategy. 

The next proposition uses the concept of a Hessian 
matrix. For any twice-differentiable function G : R" R 
and any element c = (ci, . . . , c„) € R", let 



HesSc(G) = 



dxidxj 



(c) 



(21) 



Proposition 2. Let f : {0, 1}" -> M &e a nonlocal binary 
XOR game. Suppose that the following condition holds: 
for any global maximum (ao, ai, . . . , a„) of the function 
Zf (Oq, 01, . . . ,0n), the Hessian matrix 



Hess 



(22) 



is nonsingular. Then, f satisfies a second- degree approx- 
imation property. 

Proof. Note that the Hessian condition implies that every 
maximum of the function is an isolated maximum, and 
therefore there are only a finite number of global maxima 
(afc) in the set [— tt, 7r]"+-^. 

Let s be the maximum value oi Zf (which is the opti- 
mal n-qubit score for the nonlocal game /). Additionally, 
for later use in the proof, let us define another constant. 
Consider the set 



ibnCtn 



)|} 



where we let (ao, ai, . . . , a„) vary over all the global max- 
ima for / in the set [— tt, 7r]"+^, and we let (6^) vary over 
all sequences in the set {—1,1}". The largest value in 
this set is s. Let S2 denote the second- largest value in 
this set. (If no such value exists, let S2 = 0.) 
Suppose that 



is a quantum strategy which achieves score s — e 
a unitary change of basis, we may assume that 

gift 



1 

1 











(23) 
After 

(24) 



for some (/3j). 
Let 



N= J2 /(»!,..., *«)iV} 

(U-)6{0,1}" 



(n) 



<N^-\ (25) 



We have (-0|N|?/') = s — e, and therefore ||N|| > s — e. 
Writing the operator norm of N in terms of the function 
Zf, we have 

max Z/(7,5i^i,...,&„/3„) > s-e. 

(fcfc)e{-i,i}" 
76 [— '"■.f ] 



Thus 



Zf (7, di^i, . . .,dnPn) > S - 



(26) 



for some sequence (dk) & {—1, 1}" and some 7 £ [— tt, tt]. 
By Lemma 1.3 from the Supplementary Information, 
there is a constant C (depending only on /) and a global 
maximum (ao, ai, . . . , an) for Zf such that 



\dkl3k-ak\<CV~e Vfce {l,2,...,n}. (27) 



Let 



Note that 

Let 

M = 



1 

1 



M. 



(1) _ 



M, 



(i) 



N. 



(0 



(28) 



^ /(ji,...,z„)M^)^...0M^). 

(u-)e{o,i}" 



The eigenvalues of this operator are given by 



..,6^"""") I |(6fe)e{-l,ir}, 



and its largest eigenvalue is s. 

Consider the quantity ||M — N||. Let p be the sum of 
the absolute values of the coefficients of Pf{XQ, . . . , A„). 
An easy triangle-inequality argument based on (j27p 
shows that the absolute value of any entry of M — N 
is bounded by npCy^. Therefore, 



Thus, 



IMI 



\M-N\\ <npCy/e. (29) 



> (V'|N|V')-|1M-N|1 (30) 
>s-e- npC^ft. (31) 



The quantity Si defined earlier is an upper bound for 
all the eigenvalues of M except for s. By Lemma 2.1 
from the Supplementary Information, there exists an s- 
eigenvector for M such that 



1(^10)1^ >l-(^i^^ 

S — S2 



(32) 
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and therefore 



1(^10)1 >1 



S- S2 J 



(33) 



All constants in the expression in parantheses above are 
independent of the chosen strategy (|23|) . and this expres- 
sion is O(y^) as e — > 0. The above inequality combined 
with gives the desired result. □ 

Self-testing. Let / be a binary nonlocal XOR game, 



and let (p e 
/ is a self-test for 



be an n-qubit state. Let us say that 
if the following condition holds: 



For any optimal nondegenerate n-qubit strategy 



(0) A#(l) 



My'}}j), there exist unitary matrices 
such that 



■0 = (C/i (g) U2 



..(E)Un 



(34) 



Proposition 3. Let f : {0,1}" — )■ K 6e a nonlocal bi- 
nary XOR game. Suppose that the function Zf has 
exactly two maxima in the set [—tt^tt]"''^^, and that 
these two maxima are of the form (ao, ai, . . . , a„) and 
(— QfOj ~Q!i, . • ■ , — On), with ai ^ {— 7r,0,7r} for all i G 
{1, 2, . . . , n}. Then, f is a self-test for some n-qubit state 



Proof. Let 



..0) 



P/(qi, . . . , an) 
\Pf{ai,...,an)\ 



Suppose that {-tp, {{Mj^\ M^-'^'}}j) is an optimal nonde- 
generate n-qubit strategy for /. After a unitary change 



r(i)l 



of basis, we may assume that the operators M. 
the form 



1 

1 



(1) _ 







have 



(35) 



with {9j) e [— TT, tt]", and we may make the additional 



assumption that the vectors (ai, . . . , a„) and 
lie in the same quadrant. (That is, aj > < 
for all j.) 
Again we let 



> 



M 



(ik) 



/(ii,...,i„)Mf^)® ...0M^^"). (36) 



Since the chosen strategy is optimal, by formula (pn|) we 
must have (^i, . . . , 6'„) = (ai, . . . , a„). Moreover, the 
vector must lie in the eigenspace corresponding to the 
largest eigenvalue of M. This eigenspace is spanned by 
6. The result follows. □ 



Example. Let us illustrate the previous propositions 
with the example of the Greenberger-Horne-Zeilinger 
(GHZ) paradox. Define a function h: {0, 1}^ R by 

/i(0,0,0) = 1 
/i(0,l,l) = /i(l,0,l) = /i(l,l,0) = -1 
/i(0, 0, 1) = h{0, 1, 0) = h{l, 0, 0) = /i(l, 1, 1) = 0. 

Then, 



cos 



— COS(fo 



+ 63)- COS (61c 



+ 612 + 1 



It is easily seen that Zh has two maxima inside the re- 
gion [-^,7r]4: (0,f,f,f) and (0, 



2 ' 2 ' 



-f). There- 
fore by Proposition [21 /i is a self-test for a 3-qubit state. 
(In fact, it is a self-test for the GHZ state \GHZ) — 
^ (1 000) + |111)).) Moreover, the Hessian at both of 
these maxima is equal to 



-2 -1 -1 



-2 
-1 
-1 
-2 



(37) 



which is nonsingular. 

Combining Propositions [5] and |31 we have the follow- 
ing: for any state ip G C'^ (E) <E) which is part of a 
nondegenerate strategy that achieves score 1 — e, there is 
a state ip' which is equivalent to tp under local unitary 
transformation satisfying 



\{ij',GHZ) \ > 



(38) 



where C is a fixed constant. 

The reader can compare this to the non-robust self- 
testing result of Colbeck (Q, Appendix C). 

General quantum strategies. Now suppose that we 

consider quantum strategies of arbitrary dimension. 

Whenever there are two projective binary measurements 
{p(o,-h) p(i ,-)} {p(i,+)^p(o,-)} gj^gig Hilbert 

space Q, there exists a decomposition 



(39) 



e=i 



which is respected by each of the projections P^'^*), with 
dim < 2|^ This allows us to reduce general quantum 
strategies to n-qubit strategies. Using this fact, the fol- 
lowing generalizations of Propositions [1] [21 and [31 can be 
proved. 

Proposition 4. Let f be a binary nonlocal XOR game. 
Let Q'j be the supremum of the scores achieved by quan- 
tum strategies (of any dimension) for the game f . Then, 
Q'j: is equal to the quantity Q f from Proposition [II □ 



^ See Lemma 1 from O. Regev's lecture notes, 
|wwM. cs . tau. ac . il/-od6dr/teaching/quaiituiii_f all_2005/ln/qma.pdf | 
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Proposition 5. Let f be a binary nonlocal XOR game 
which satisfies the condition of Proposition Then, 
there exists a constant C > such that the following 
holds. For any quantum strategy 

ip e Qi <E> . . . (g> Qn, (40) 
{Mf : Q,^Q,}^^. (41) 

whose score is within e of the optimal score s, there exists 
a quantum strategy 

(/)e7^l(g)...0 7^„, (42) 

{iV«:7^,^7^,}^^ (43) 

achieving score s and unitary embeddings 

Uf.Q.^n, (44) 

such that 

\{{Ui® ...(^Un)^,^)\>l-C^e. □ (45) 

Proposition 6. Let f be a binary nonlocal XOR game 
which satisfies the conditions of Proposition O Then, 
there exists an n-qubit state a € (C^)*^" such that the 
following holds: for any n-partite state ip d Qi (8) . . . ® Qn 
which is part of an optimal strategy for f , there exists 
another n-partite state X € Qi . . . ^.''^d unitary 

embeddings Vj : Qj ^ ® Q'j such that 

(Fi (g) . . . (g) Vn) IV') = |a) ® |A) . □ (46) 

Conclusion. We have provided a general scheme for 
proving robust self-testing results in the context of binary 
XOR games, and we have applied it prove a new result 



on the GHZ paradox. A natural question is whether the 
scheme can be generalized beyond binary games. Such 
generalization may be difficult since some of the mathe- 
matical tools we have used (such as the canonical form 
([T5|) and the decomposition ([Ml)) break down when non- 
binary inputs and outputs are allowed. Some new math- 
ematical insight may be required to extend our results. 

It would be interesting to try to apply the GHZ re- 
sult to randomness expansion, following proposals such 
as those in Q and 
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SUPPLEMENTARY INFORMATION 



1. Lemmas for multivariable functions 



For any twice-differentiable function F : 
let HesSci^ denote the Hessian matrix, 



and any c = (ci, . . . , c„) G 



HesScf := 



Lemma 1.1. LetG: 



ozf 
dz2 dzi 



-F 



^ 

dzi dz2 

Kf 



dZn dz2 ^ 



dzi dzn 

02:2 dZn 



rF 



(ci,...,c„). 



— !■ K 6e a twice-differentiable function. Letn. — (xi, . . . ,a;„) G 
R" he such that all first order partial derivatives of G vanish at x, and the Hessian 
matrix HesSx(G) is negative definite. Then, there exist constants (5o,C > such 
that for all y G 

(1) 



with ||y — x|| < So, 

G(y) <G(x)-C||x-y|| 



Proof. Let A be the eigenvalue of HesSx(G) that is closest to zero, and let C = —A/2. 
Define a function : R" ^ R by 

(2) ^y) = G(y) + G||y-x|r 
Then, the Hessian of at x is 

(3) HesSx(G) + GI„, 

which is negative definite. Therefore, x is a local maximum of F. The result 
follows. □ 

For any vector y = (yi, . . . ,y„) G M", let ||y||^ = max^ \yi\. Note that 

(4) l|y|loo< llyll < V^llylloo- 

The following modification of Lemma 11.11 follows immediately. 

Lemma 1.2. LetG: R" — > R 6e a twice-differentiable function. Letx. = {xi, . . . ,Xn) G 
R" be such that all first order partial derivatives of G vanish at x, and the Hessian 
matrix HesSx(G) is negative definite. Then, there exist constants 5o,G > such 
that for all y G R" with ||y — x||^ < Sq, 

(5) G(y) <G(x)-G||x-y||^. □ 
For any set S" G R" and any (5 > 0, let 

(6) BiiS) = {y G M" I 3x G 5 such that |lx -y\\^<S}. 
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Lemma 1.3. Let G : M" — > M &e a twice- differentiable function. Let c be a positive 
real number. Let b be the maximum value of G achieved on the set [— c, c]". Suppose 
that this maximum is achieved only by a finite number of points jx^^-', . . .x*-™-*}, 
all of which occur in the interior of [— c, c]". Suppose also that 

(7) detHess^(.)(G) 7^ 

for all i e {1, . . . , m}. Then, there exists a constant K > such that for all e > 0, 

(8) G-i([6-e,6])CS^V^(y). 



Proof. By Lemma 11.21 there exist constants C, eo > such that for any i G 
{1, . . . , m} and any y G M" such that lly — x^'-' 11 < cq, 



I oo 

2 



(9) G(y)<G(x«)-C xW-y 
Let b' be the largest value achieved by G on the compact set 

(10) (ye[-c,c]"| y-x« >eo Vz £ {1, . . . , n}| . 
Setting 

(11) K = nmxlc-\^L=\ 

yields the result. □ 

2. A LEMMA FOR HeRMITIAN OPERATORS 

Lemma 2.1. Lei X is an n x n Hermitian operator. Suppose that r is the largest 
eigenvalue of X , and suppose that r' is a quantity that is less than r but greater 
than or equal to all other eigenvalues of X . Suppose that v is a unit vector such 
that 

(12) {v\X\v)>r-e. 

Then, there exists a unit r-eigenvector u for X such that 

(13) |(^|„)|2>i_^. 

r — r 

Proof. Let V denote the r-eigenspace for X, and let denote its complement. 
Let 

(14) X' ^r ■ projv + r' ■ projy_L . 

Then, X' >X, and thus 

(15) {v\X'\v)>{v\X\v)>r-e. 
Since 

(16) {v\X'\v) - rWpvoiyvW^ +r'\\^voiy^v\\^ 

(17) = r ||projyu||^ +r' (^1 - ||proj^^i>| 
we therefore have 



(18) r llprojv^i-f + / (l - \\vYoiyv\\- 
By an algebraic manipulation, 

(19) ||proj,,i;f > 1 - ^ 



> r — e. 
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Letting u — projyw/ ||projyz;|| yields the result. □ 



Randomness Expansion from the 
Greenberger-Horne-Zeilinger Paradox"* 

Carl A. Miller Yaoyun Shi 



A Randomness Expansion Devices 

(~| ■ In this section, we formalize the notion of a multi-part randomness expansion 

PlI (RE) device. 

A 3-part randomness expansion (RE) device D is a device consisting 
of 3 components, Di, D2, -D3, satisfying all of the following conditions: 

1. Each component Dj accepts a single bit ij e {0, 1} as input and returns 
a single bit Oj € {0, 1} as output. 

2. Each component Dj contains a quantum system, Qj. When it receives its 
input bit ij, it performs an orthogonal binary measurement on Qj. The 
output bit Oj is the result of this measurement. 

3. The choice of measurement performed by Dj is determined solely by its 
input bit ij. In particular, it does not depend on the inputs or outputs of 
any of the other components. (The components do not communicate with 
one another.) 

We will say that the input bits (ij) and the output bits (oj) pass the GHZ 
test if the following Boolean relations are satisfied: 

H®i2®»3 = 0, (1) 
Oi®02©03® (Zi Vi2V«3) = 1. (2) 

In what follows, we will often be choosing the inputs to device D according 
to a random probability distribution. In this case, we denote the inputs by the 
random variables /i, I2, and /3, and the outputs by the random variables Oi, 



*This appendix was written during June - August 2011. The text given here is unmodified 
except for a few minor edits. 



1 



02, and 03- Let / = {h,h,h), and O = (Oi, 02,03). 

h h h 





Qi 








Q2 




D2 




Q3 





D3 



Oi 



02 



03 



Let Q denote the tripartite quantum system formed by Qi, Q2, and Q3. As a 
convention, if Z is a quantum system, then we use the corresponding symbol Z 
to denote the complex Euclidean state-space of Z. Thus Qi, Q2, and Q3 denote 
the state spaces of Qi, Q2, and Q3. 

Note that the state of the system Q may be mixed. Sometimes we will 
consider an additional quantum system that is entangled with Q: 

h h h 



D2 



D3 



E 



Oi 



O2 



O3 



We will use the expressions rl^'^'^^ and rl^°**^ to denote the pre- and post- 
measurements states of the systems above. Thus, 



-,(post) 
IQE 



(3) 



denotes the post-measurement state of the systems 7, Q, and E, taken together. 
We will drop the "(pre)" and "(post)" superscripts when they are not necessary. 

r is always presumed to denote a density operator. If the state of a collec- 
tion of quantum systems happens to be pure, then we may also also use the 
symbol v to denote the vector representing the pure state. (Thus for example, 
if the expression Uq"^^^ appears, it should be interpreted to mean a vector in Q 
representing the pure prc-measurement state of Q.) 

A.l A canonical form 

Let D he a. 3-part RE device (as defined above). For any j G {1, 2, 3}, and any 
input bit i G {0, 1}, let 



(4) 



2 



be the orthogonal measurement the the component Dj performs on its quantum 
system Qj on input i. This is a set consisting of two complementary orthogonal 
projection operators on the space Qj. 

It is convenient to express as a single operator. Let 

^ pirn _ pi^A) ^ (5) 

This is a Hermitian operator on Qj whose square is equal to I. In general, any 
Hermitian operator X satisfying X'^ — 1 yields a binary orthogonal measurement 
({i^, ^--^})- For this reason, we will often refer to such an operator X simply 
as a "measurement." 

Let D' be another 3-part RE device. 



Q'2 



whose measurements are denoted by T'^j K Let us say that a unitary embed- 
ding of D into D' is a collection of unitary embeddings 



*3 



Qi S'l 

Q2 Q'2 
Q3 S3 



which map the pre-measurement state of Q = (Qii Q2, Qs) to the pre-measurement 
state of Q' = (Qi, Q's): ^-nd which are such that 



(6) 



It is easy to see that if such a unitary embedding exists, then the device D' has 
the same output-statistics as D. 

If E is an additional quantum system, and the pre-measurement states of D 
and D' are given as entangled states 



QE 



and F'^^; 



(pre) 



E ' 



then we say that the unitary embedding {^i) respects entanglement with 
E if the superoperator 



(*1 (g) «) *3 (g) Is) (•) (*i (g) ^'2 «> ^'3 «> If) 



(7) 



carries F 



(pre) 
QE 



to F' 



(pre 



Proposition A.l. Let D he a S-part RE device which is entangled with an 
additional quantum system E. Then there exists a unitary embedding of D into 
another 3-part RE device D' , respecting entanglement with E, such that the 
following conditions hold. 



3 



1. The state-spaces of D' have the form 

Q'l = A® 111 
Q'2 = B® 7^2 
Q'^ = C(8)7^3. 



(8) 
(9) 
(10) 



where A, B, and C are complex Euclidean spaces, and TZi, TZ2, o,nd TZs 
are copies of C^^'^^ . 



2. The measurements operators for Dj on input i — are given, respectively, 
by 

' 1 

1 



1 

1 

1 

1 



^ |afe) (flfel ® 

k 

E (^^1 ® 

E 1"^™) ® 

m 

where {ok}, {be}, and {cm} denote orthogonal bases for A, B, and 
3. The measurements operators for Dj on input i — 1 are given by 

E l"*^^ ^'^'=1 

k 



(11) 
(12) 
(13) 



'-■m/ \<--m I 






Afe 










7i 


_ le 








</>m 


<j)rn 






(14) 
(15) 
(16) 



where Xk, 7£, and (pm denote complex numbers that satisfy 

|Afc| = \je\ = |(/)„i| = 1, 
Im(Afe) > 0, 
Ini(7£) > 0, 
lni(0„) > 0. 



(17) 
(18) 
(19) 
(20) 



Proof. It is easy to show that, given any two Hermitian operators X and Y on 
such that X'^ = Y^ = I, there is a unitary transformation J7 : — such 
that 



uxu"^ = 

UYU'' = 



1 

1 

C 



(21) 
(22) 



^The 2x2 matrices here are written with respect to the computational basis {|0) , |1}} of 



4 



where ^ satisfies |C| = 1 and Im(C) > 0. From this fact, it is easy to prove 
the proposition for the case where the spaces Qj aU have dimension < 2. The 
general case then fohows via Jordan's lemma0 □ 



B The GHZ Paradox: 2x2x2 case 

We are interested in those RE devices which pass the GHZ test (see ([2])) with 
high probabihty. Because of Proposition lA.li it is useful to focus on RE devices 
that are in a particular form. 
Let J be a 3-part RE device. 



Ri 



Ji 



i?2 



J2 



J3 



Suppose that TZj = C^"'^^ for each j, and suppose that the measurement oper- 
ators {5* -'^} for J are given by 



c(0) _ c(0) _ oiO) 
'-'1 — '-'2 — '-'3 



(1) 



1 

1 

A 

A 

7 

7 

(j) 

6 



(23) 
(24) 
(25) 
(26) 



with |A| = I7I = \(/)\ — 1. Suppose that the pre-measurement state of (i?i, R2, R3) 
is a pure state given by a vector a e TZi TZ2 ^ TZ^. 

Consider that a 3-bit string /1/2/3 is chosen uniformly at random from the 

set 



{000,011,101,110}, 



(27) 



and given to device J. Then J passes the GHZ test if and only if the output 
string O1O2O3 satisfies the relation 



Oi © O2 ® O3 © (/i V /2 V I3) = 1. 



(28) 



We can calculate the probability that this formula will be satisfied using the 
operators Sj^\ For example, the probability that the formula will be satisfied 
when the input string is 110 is: 



1 



2^i^i 



(1) 



,(1) 



:.(o) 



) a, a) 



(29) 



^See Lemma 1 from |http: //mmw. cs ■ tau. ac ■ ll/-odedr/teaching/quantum_f all_2005/lii/qma.pdf | 



The probabihty that this formula will be satisfied on an input randomly chosen 
from (Ull) is 



a 



(30) 



Let 



a= ^ Ckim \klm) 

k,l,me{0,l} 



(31) 



Writing (1301) in terms of A, 7, 0, and {cjfci}, and rearranging terms, we find the 
following expression: 



P (J passes the GHZ test) 



Re 
Re 
Re 
Re 



Cm ■ cooo 



Clio • cooi 



cioi ■ coio 



ClOO • Coil 



-1 - 


h A7 + 70 H 


^0A 




4 




-1 - 


h A7 + 70 H 


-0A 




4 




-1 - 


h A7 + 70 H 


h0A 




4 




-1 - 


h A7 + 70 H 


^0A 



We will make frequent use of this expression in the subsections that follow. 
B.l A Few Lemmas 

The following lemmas will be useful in subsequent sections. 

Lemma B.l. Let a, h, and c be complex numbers such that \a\ — \b\ = \c\ — 1 
and Im(a) > 0, lm{b) > 0, and Im(c) > 0. Then, 



'1 + at + be + ca 



< \ I 



(32) 



Proof. Consider, for all complex numbers b' satisfying = 1, the value 

\-l + ab' + b'c + ca\ = + ca) + b' {a + c)\ . (33) 

This value will be maximized when the angle of (— 1 + ca) in the complex plane 
agrees with the angle of b'{a + c). This occurs precisely when b' = i. Therefore, 



\ — l + ab + bc + ca\ < \ — l + ai + ic + ca\ 



(34) 
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By similar reasoning, 

| — 1 + ai + ic + ca\ < I— 1 + ai + + ml . (35) 



Therefore, 



-l + ab + bc + ca\ < |-2 + 2ai\ = \2a + 2i\ (36) 



Note that \2a + 2if + \2a - 2if = 16. Therefore, 



imphes 



\-l + ab + be + caf < \2a + 2if (37) 



-1 + ab + be + ca\^ < 16-\2a-2if. (38) 



Dividing by 16 and taking square roots of both sides of the equation yields the 
result. □ 

Lemma B.2. Let a, b, and e be eomplex numbers sueh that \a\ = \b\ = \e\ = 1. 
Suppose that Im(a) > and Im(6) > 0, but Im(c) < 0. Then, 



-1 + ab + be + ca 



< f • (39) 



4 

Proof. Consider the complex numbers 

i-l + ab) and (a + 6). (40) 

In the complex plane, (— l + a6) lies at an angle of +7r/2 (in the counterclockwise 
direction) from [a + b). Since Im(c) < 0, the angle between the product c{a + b) 
and (— 1 + ab) must be an obtuse angle or a right angle. Therefore, 

\{-l+ab) + c{a + b)\^ < |-1 +a6|^ + \a + bf < 8. (41) 

The lemma follows. □ 

Lemma B.3. Suppose that A, B, A' and B' are Hermitian operators on C" 
whieh satisfy \\-\\^ < 1. Then, 

\\A^B-A'^B'\\^ < \\A-A'\\^ + \\B-B'\\^. (42) 

Proof. We have the following: 

\\A<g>B-A'<^B'\\^ = \\(A-A')^B + A'^{B-B')\\^ (43) 

< ll^-^'llooll^lloo + P'ILII^-S'lloo (44) 

< \\A-A'\\^ + \\B-B'\\^. (45) 

□ 
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B.2 Characterizing Approximate GHZ Devices 

There exists 3-part RE devices which satisfy the GHZ test perfectly. Consider 
the 3-qubit device whose pre-measurement state is the pure state 



g=^ (|000) - |111)) e C^°'^> ® C{°'i> C{°'i> 



(46) 



and which is the such that each component uses the measurement operator 



1 

1 



(47) 



on input i = 0, and the measurement operator 



i 
-i 



(48) 



on input i = 1. This device passes the GHZ test with probabihty 1. Let us refer 
to this device as the ideal GHZ device. 

We will use the above example for comparison. In the next two propositions, 
we see that if a 2 x 2 x 2 RE device passes the GHZ test with high probability, 
then it is similar to an ideal GHZ device. 

Proposition B.4. Let J be a 3-part RE device whose state spaces TZj are copies 
of the qubit-space C^^'^^ , and whose pre-measurement state is a pure state 



a= ^ Ckim \klm) . 

k,l,me{0,l} 



(49) 



Suppose that the measurement operators for J are 



(0) 



5[ 



1 

1 

A 

A 



5f = 



1 

1 

7 

7 



4°^ = 



^>3 — 



1 

1 

6 



(50) 



where 



|A| = |7| = - 1 
Im(A) > 
Im(7) > 
Im((^) > 0. 



(51) 
(52) 
(53) 
(54) 



Suppose that 



P (J passes the GHZ test) = 1 — e. 



(55) 
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Then, 



\X-if < 16e 
I7 — if < 16e 
\(b-if < 16e. 



(56) 
(57) 
(58) 



Proof. From the discussion at the beginning of this section, we know that the 
probabiHty that J passes the GHZ test is given by 



+ Re 

+ Re 

+ Re 

+ Re 



Cm • cooo 



Clio • cooi 
cioi • CoiO 
ClOO • Coil 



-l + A7 + 7(/) + (/)A \ 
-l + A7 + 7(?^ + ^A \ 

4 _ _ ; 

—l + Xj + j(j) + (j)X\ 



Therefore, the quantity (1 — e) satisfies 



1-e < 2 + |ciii||cooo| 



+ |ciio||cooi| 
+ |cioi||coio| 
+ |cioo||coii| 



-1 + X-f + ^(j) + (f>X 



-1 + Xj + -f(p + (j)X 



-l + A7 + 70 + (?iA 



- 1 + A7 + 7^ + (/)A 



which imphes 

1-e < 



+ 



ciiil^ + |cooo|^) 
|cooi|^) 



-l + A7 + 7^ + ^A 



+ (|ciio|V 
+ (|cioi|^ + |coio|^) 



1 + X'j + ^(j) + (j)X 



-1 + A7 + 7(/) + ^A 



+ 



(|cioo|^ + |coii|^) 



- 1 + A7 + 7<?!) + (?!>A 



(59) 



(60) 



(61) 
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or equivalently, 

1 - 2e < 



(|cill|^ + |cooo|^ 



ciiol 



ClOll 



Clool 



|C00l| 



|coio| 



|coii| 



-1 - 


1- A7 + 


7(/)H 


L0A 




4 






-1 - 


h A7 + 


7^H 


^0A 




4 






-1 - 


F A7 + 


7(/)H 


^<j)X 




4 






-1 - 


F A7 + 


7^H 


^0A 


4 



(62) 



Consider the following four quantities from inequality (|62 



-1 + A7 + 7(/) + 0A 



-1 + A7 + 7(/) + (/)A 



-1 + A7 + 70 + 0A 



-l + A^ + T^^ + cZ-A 



Note that the coefficients of these quantities in ((62)) sum to 1. By Lemma [6.21 
the second, third, and fourth quantities are each less than or equal to \/2/2. If 
the first of the quantities were also less than or equal to a/2/2, then we have 
1 - 2e < and therefore e > ^ - ^ > |, and inequalities ([55]) - dSH) would be 
trivially satisfied. So, let us assume that the first of these quantities is greater 
than a/2/2, which makes it the largest of the four. Then, 



1 - 2e < 



-1 + Xj + + (j)X 



By Lemma FB. II 



Therefore, 



1 - 2e < 



(1 - 2ey < 1 



\X- 



1 _ 4g < 1 _ 



iA-*r 



16e > \X-if 

Similar reasoning shows that I7 — < 16e and \(f> — < 16e. 

Proposition B.5. Let J be the device from Proposition \B.4\ Let 

1 



Then, 



- ^(1000)^1111)) 
|(a,5)|>l-7e. 



(63) 

(64) 

(65) 

(66) 
(67) 
□ 

(68) 

(69) 
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Proof. As before, we know that the probabihty that the probabihty that J 
passes the GHZ test is given by (|59p . We prove the proposition in 3 cases. 



Case 1: Cqoo = £111 = ^■ 

Using Lemma IB.21 we have the fohowing: 



Therefore, 



< 



< 



Clio • cooi 



ClOl • CoiO 



Cloo • Coil 



Re 
Re 
Re 
|ciio| |cooi| 
|cioi| |coio| 
|cioo| |coii| 



|ciio|^ + |cooi|^ 



-1 - 


h A7 + 


70 H 


-0A 




4 






-1 - 


h A7 + 


70 H 


-0A 




4 






-1 - 


h A7 + 


70 H 


-0A 



-1 + A7 + 70 + 0A 



4 






-1 + A7 + 


70 H 


^0A 


4 






-1 + A7 + 


70 H 


h0A 





2 




|cioi|^ 


+ 


CoiO^ 




2 




|cioo|^ 


+ 


|coii|^ 


2 



4 

V2 
2 

V2 
2 

V2 
2 



V2 
4 



e > 1 



(70) 



(71) 



(72) 



(73) 



Since the quantity above is greater than y, condition ([M)) is triviahy satisfied. 
Case 2: Cqqi = Cqiq = Cioo = Cqh = Cioi = Cuo = 0. 
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We have 



1 - e = ^ + 



Cm • cooo 



= Re 

< Re 

< Re 

< Re 



Cm • Cooo 



Cm • Cooo 



Cm • Cooo 



Cm • Cooo 



-1 

T 

-1 
T 

-1 
T 

-1 
T 



Cm • Cooo 



+ Re 

+ Iciiil Icooo 



I 4 



A7 + 7(/) + (/)A 



+ 



|ciii|^ + |cooo|^ 



which imphes, by a hnear manipulation, 

1 - 8e < (-2) • Re [cTTTcooo] 

We have the following: 

2 - 8e < 1 - 2 • Re [cTTTcooo] 



(74) 

(75) 

2 



By an easy calcuation, the right side of this inequality is equal to 2 |(a,5)| 
Thus, 



2- 


8e 


< 


2|(«,3)l' 


1 - 


4e 


< 




1 - 


4e 


< 




1 - 


7e 


< 





(76) 
(77) 
(78) 
(79) 



Case 3: General case. 

Let 

a = xa' + ya" (80) 

where a', a" are unit vectors satisfying the conditions of Case 1 and Case 2, 
respectively, and {x,y) is a unit- length vector in R^. Let Jq,/ and Ja" denote 
the device J with the state a replaced by a' and a", respectively. It is clear 
from expression (|59|) that 

P (J passes the GHZ test) — x'^ -P {Ja' passes the GHZ test) 

+ • P ( Jq" passes the GHZ test) . 

This case follows from cases 1 and 2 by an easy linearity argument. □ 
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Corollary B.6. Let J be the device from Proposition \B.4\ There exists a com- 
plex number ( with |C| = 1 such that 

||a-C5ll'<14e. (81) 

Proof. Choose C so that the inner product {a, Qg) is a nonnegative real number. 
(Thus, (a,C5) = \ {a,g)\.) Then, 



l|a-C.9ll = {a,a) - {aXg) - {C9,a) + (CgXg) 

= 2-2|(a,5)| 

< 2-2(l-7e) 

= 14e. 



(82) 
(83) 
(84) 
(85) 

□ 



B.3 The Post-Measurement State of an Approximate GHZ 
Device 



Let J be the 3-part RE device from Proposition IB.4I Suppose, as usual, 
that J is given an input string /1/2/3 that is randomly chosen from the set 
{000,011,101,110}. Consider the post-measurement state of J together with 
its input and output registers. 

h h h 



Ri 



Ji 



R2 



J3 



Ox 



O2 



O3 



Recall that the pre-measurement state of i? = (i?i, i?2, -R3) is a pure state a. 
We can express the post-measurement state of lOR as a pure state like so: 



(post) 
"lOR 



I ) ' X! I*l«2«3) (8) I01O2O3) ( 



(86) 



(ifc),(ofc) 



'l + (-l)°iS'J'^^\ /l+(-l)°^4 



where the summation is taken over all input /output combinations (ik), (o/c) 
which satisfy the GHZ test (([I])-©). 
Let 



5--^(|000)~|lll)) 



(87) 
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and 



1 

1 



.(1) 



i 
-i 



(88) 



Define a new state, denoted ^^/ofl'''' expression ([551) with the operators s\^^ 
replaced by a*^*) and the state a replaced by g: 

(ideal) _ 



"lOR 



• Nl*2«3) «) |0l0203)(8) (89) 

(ifc),(ofc) 

I + (-l)°icr(*i) \ /I + (-l)°2cr(*i' \ /I + (-l)°-^cr('i) ■ 

■ 5- 



This is the post-measurement state of the ideal GHZ device from subsection lB.21 
Let us compare expressions (1551) and (15^ . From Proposition IB. 41 we have 



- a^^-^) 



< 32e. 



(90) 



Therefore, 

oo 2 

< \/32i. (92) 
The individual projection operators from (1551) and (15^]) therefore satisfy 



'l + (-l)'"=S'f \ /I + (-l)°'=cr(»'=) 



< 



/32e 



(93) 



Let M and m(*'^'="') be the operators enclosed in brackets in ^ and (|5g)) . 
respectively. Applying Lemma IB. 31 we find that 



< 3 • V8e. 



By Corollarv lB.61 there is a unit-length complex number C, such that 



(94) 



(95) 



Therefore: 



(post) 
■'lOR 



(ideal) 
lOR 



< 
< 

< 3 



Ma- M'^"^'"'^'>{Cg) 

Ma- M^'^^'^'^'^a + m'^"^""^^ a - M'^"^'"'^\(:g) 



+ l|M|L||a-C5|| 
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Squaring, we find 



(post) 



'lOR 



> (ideal) 



2 

14) e 



(96) 
(97) 



< 

< 150e. 

We state the above inequality as a proposition. 
Proposition B.7. Let J be the device from Proposition \B.4\ Suppose that 

P (J passes the GHZ test) = 1 - e. (98) 



Then, there exists a unit-length complex number ( such that 

2 



(post) 
lOR 



J. (ideal) 
' ^lOR 



< 150e 



(99) 



where d"''' denotes the post-measurement state of the ideal GHZ device. 



C The GHZ Paradox: General Case 

Now we will broaden our focus and consider 3-part RE devices that contain 
arbitrarily large quantum systems. 



C.l The Post-Measurement State of an Approximate GHZ 
Device: General Case 



As in subsection lB.31 if i? is a 3-qubit system, then we write 



(ideal) 



for the post- 



measurement state of an ideal GHZ device on R. Also, we will write F^q^'^ for 
the corresponding density operator. 

Theorem C.l. Let D' be a 3-part RE device that is in the canonical form 
described in Proposition \A.1[ Let F be a quantum system which is entangled 
with D' in such a way that the state Tqi p is pure. Let L and O denote the input 
and output registers of D' . 



h 



h 



h 



Q'2 



D' 



F 



D' 



Suppose that 



O2 



O3 



P {D' satisfies the GHZ paradox) = 1 — e. 



(100) 
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Then, there exists a unit vector 

w e A (g) B -g) C (g> T 

such that the inequality 



(post) 



(ideal) 



< 150e 



(101) 



(102) 



holdsE 



Proof. Choose an orthonormal basis {/„} for F. Let us write the pre-measurement 



state Vq^p as 



^Q''f = 51 (a/c ® ^/ c„ (K) /„) (8) akir, 



(103) 



klmn 



with aklmn G 7?.i (g) 7?.2 (g) 7?.3 . 

For every 4-tuple {k,l,m,n) which is such that akimn 7^ 0, define Jkimn to 
be the 3-part RE device on the system {Ri, R2, R3) whose pre-measurement 
state is 



I ^klmn I 

and whose measurement operators are given by 



(104) 



SI 



(1) 






1 


1 








Xk 








5f = 



1 

1 

7/ 

77 



c(0) _ 
— 



S. 



(1) _ 



1 

1 

0m 





(105) 



For each device Jkimm choose a unit-length complex number Ckimn so as to min- 
imize the Euchdean distance between Ckimn • ^'io'r^'^ ^^'^ ^^"^ post-measurement 
state of Jkimn- Let 



w 



^ (afc ®bl®Cm® fn) <^ C 



kl ran K^fc^mn 



(106) 



klmn 



Note that 

P {D' passes the GHZ test) — ^ |afc;mn|^ P {Jkimn passes the GHZ test) 



klra\ 

Likewise, we can express the quantity 



{ideal) 
lOR 



(107) 



•^Note that in expression II102I I. we have treated the vector spaces TZj from the statement 
of Proposition I A. II as the state spaces of a 3-qubit system. The expression 'Wjq^'' represents 
the post-measurement state of an ideal GHZ device on this system. 
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as a weighted sum of the squares of the distances between the vectors CfcZmn ■ v 



{ideal) 
lOR 



and the post-measurement states of Jkin 
Proposition IB. 71 by hnearity. 



Inequahty ()102p thus fohows from 

□ 



Let us consider how Theorem lC.ll mav be strengthened. We wish to cover the 
case in which the stated device D' is entangled with another quantum system Fi , 
but the state ^'ql'p^ is not pure. In this case, we can simply introduce a second 

entangled quantum system F2 such that ^'q^p^p^ is pure. By the theorem, we 
find that for some vector w in the tensor product of A, B, C, J-i, and T2, 



(post) 
IOQ'FiF2 



W ® V 



(ideal) 
lOR 



< 150e. 



(108) 



We would like a similar expression with the system F2 omitted. Since we cannot 
use vector-length to measure distance in this case, we use the trace-norm instead. 



lOQ'Fi 



-^{ideal) 
lOR 



< 

< 
< 



[post] 
IOQ'FiF2 



lOR 



^ (ideal) 
W (g) v\oR 



600e. 



We have the following alternate statement of the theorem. 

Theorem C.2. Let D' he a device which is in the canonical form described 
in Proposition \A.1\ Let Fi be a quantum system which may be entangled ( in 
a possibly mixed state) with D' . Suppose that D' passes the GHZ test with 
probability 1 — e. Then, there exists a density operator ^abcFi on A(^B(^C(^J-i 
such that 



(post) 
lOQ'Fi 



(ideal) 
lOR. 



< 600e. 



(109) 



Secondly, we would like to have a version of the theorem which applies to 
devices that are not in canonical form. To simplify matters, let us drop the 
quantum system Q' and all of its subsystems from ()109|) . We have: 



^(post) 
lOFi 



Fi 



(ideal) 
10 



< 600e 



(110) 



Since a unitary embedding of one entangled 3-part RE device into another 
obviously has no effect on the quantity in this inequality, we have the following 
by Proposition lA.il 

Theorem C.3. Let D be a 3-part RE device (not necessarily in canonical form) 
which passes the GHZ test with probability 1 — e. Let E be a quantum system 
which may be entangled with D. Then, there exists a density operator ^p on £ 
such that 



^(post) 
lOE 



^(ideal) 
10 



< 600e. 



(Ill) 
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Finally, note that inequality (jllip above implies 



\Te - ^bH < 600e. 



Therefore, 



^{ideal) 
lO 



^{ideal) 
lO 



< 600e. 



Applying the triangle inequality for 1 1 • 1 1 , we find that 

2 



(post) 
lOE 



lO 



< 2400e. 



We therefore have the following final version of our theorem. 



(112) 



(113) 



(114) 



Theorem C.4. Let D be a 3-part RE device which passes the GHZ test with 
probability 1 — e. Let E be a quantum system which may be entangled with D. 
Then, 



-^(post) 
lOE 



^{ideal) 
lO 



< 2400e. 



(115) 
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